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A general theory for the existence of solitary wave and double layer at M — M c has 
been discussed, where M c is the lower bound of the Mach number M, i.e., solitary 
wave and/or double layer solutions of the well-known energy integral start to exist 
for M > M c . Ten important theorems have been proved to confirm the existence of 
solitary wave and double layer at M = M c . If V(0)(= V(M, 0)) denotes the Sagdeev 
potential with is the perturbed field or perturbed dependent variable associated 
with the specific problem, V(M, <fi) is well defined as a real number for all M G M. 
and for all G $, and V(M,0) = V'(M, 0) = V"(M c ,0) = 0, V"'(M c ,0) < 
(V"'(M c ,0) > 0), OV/dM < for all M(e M) > and for all 0(g $) > (0(g 
$) < 0), where " ' = d/d<fi ", the main analytical results for the existence of solitary 
wave and double layer solution of the energy integral at M — M c are as follows. 
Result-1: If there exists at least one value Mq of M such that the system supports 
positive (negative) potential solitary waves for all M c < M < M , then there exist 
either a positive (negative) potential solitary wave or a positive (negative) potential 
double layer at M — M c . Result-2: If the system supports only negative (positive) 
potential solitary waves for M > M c , then there does not exist positive (negative) 
potential solitary wave at M — M c . Result-3: It is not possible to have coexistence 
of both positive and negative potential solitary structures (including double layers) 
at M — M c . Apart from the conditions of Result-1, the double layer solution at 
M = M c is possible only when there exists a double layer solution in any right 
neighborhood of M c . Finally, these analytical results have been applied to a specific 
problem on dust acoustic waves in nonthermal plasma in search of new results. 
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I. INTRODUCTION 



The Sagdeev potential approach^ has been considered to investigate the existence of 
solitary wave and double layer at M — M c , where M c is the lower bound of the Mach 
number M, i.e., solitary wave and/or double layer solutions of the energy integral start 
to exist for M > M c . In most of the earlier works,- - — solitary wave and/or double layer 
solutions have been investigated for M > M c . However, some recent investigations^"— have 
shown that finite amplitude solitary wave can exist at M = M c in the parameter regime 
where solitons of both polarities exist. The numerical observations^"— of the solitary wave 
solution of the energy integral at M = M c , motivate us to set a general analytical theory for 
the existence of solitary wave and double layer solution of the energy integral at M — M c . 
The present paper is purely analytical study on the existence of solitary wave and double 
layer solution of the energy integral at M — M c . Without this analytical study, it is difficult 
to predict the existence of double layer solution at M — M c . The analytical study also give 
the point in the parameter regime where the double layer solution at M = M c exists. So, we 
consider the following well known energy integral, which has been commonly used in various 
plasma environments.-"- 1 ^"— 



with V{4>) = V(M,<p) is the Sagdeev potential. Here, V(M, 0) is same as V(<fi), i.e., the 
Mach number M is omitted from the notation V(M, 0) when no particular emphasis is put 
upon it and let M. x $ is the domain of definition of the function V(M, 0), i.e., V(M, 4>) is 
well defined as a real number for all M G M. and for all <p e Hereafter, by the phrases " 
for all > " , " for all < " , and " for all M > " , we mean, respectively, " all those 
G $ which are strictly positive " , " all those G $ which are strictly negative " , and " all 
those MeM which are strictly positive " . 

If " ' " indicates differentiation of V(M, 0) with respect to 0, then with the help of some 
definitions and some well-known theorems of real analysis, we have proved the following 
two theorems, which are necessary to prove the main theorems regarding the existence of 
solitary wave and double layer at M — M c . 



Theorem 1 : If V"'(M C , 0) ^ 0, then there exists a strictly positive real number e (> 0) 
such that V'"(M C , 0)V m (M e , 0) > for all -0 e < < e . 




(1) 
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Theorem 2 : If V"(M C , 0) = and V'"(M C , 0) ^ 0, then there exists a strictly posi- 
tive real number e (> 0) such that V"'(M C , 0)V"(M c , 0) > for all < < e and 
V"'(M C , 0)V"(M c , 0) < for all -<f> e < < 0. 

If V(M,0) = V'(M,0) = V"(M c ,0) = 0, using Theorem 1, Theorem 2, and some 
well-known theorems of real analysis, we have proved the following theorems to confirm the 
existence of solitary wave and double layer solution of the energy integral (fTT) at M — M c . 
Theorem 3 : If V"'(M C , 0) < 0, then at M — M c , there may exist either a solitary wave 
or a double layer in the positive potential side but there does not exist any solitary wave or 
double layer in the negative potential side. 

Theorem 4 : If V"'(M C , 0) > 0, then at M — M c , there may exist either a solitary wave or 
a double layer in the negative potential side but there does not exist any solitary wave or 
double layer in the positive potential side. 

Theorem 5 : If V"'(M C , 0) ^ 0, then it is not possible to have coexistence of both positive 
and negative potential solitary structures (including double layers) at M = M c . 
Theorem 6 : If V'"{M C , 0) < 0, dV/dM < for all M > and for all > 0, and if 
there exists at least one value M of M such that the system supports Positive Potential 
Solitary Waves (PPSWs) for all M c < M < M , then there exist either a PPSW or a Positive 
Potential Double Layer (PPDL) at M = M c . 

Theorem 7 : If V"'(M C , 0) > 0, dV/dM < for all M > and for all < 0, and if there 
exists at least one value Mo of M such that the system supports Negative Potential Solitary 
Waves (NPSWs) for all M c < M < M , then there exist either a NPSW or a Negative 
Potential Double Layer (NPDL) at M = M c . 

Theorem 8 : If V'"(M C , 0) < 0, dV/dM < for all M > and for all > 0, and if 
the system supports only NPSWs for M > M c , then there does not exist any PPSW at 
M = M c . 

Theorem 9 : If V'"(M C , 0) > 0, dV/dM < for all M > and for all < 0, and if 
the system supports only PPSWs for M > M c , then there does not exist any NPSW at 
M = M c . 

Theorem 10 : If V"'{M C , 0) ^ 0, dV/dM < for all M > and for all 0^0, and if 
there exists at least one value Mo of M such that the system supports both PPSWs and 
NPSWs for all M c < M < M , i.e., if the system supports coexistence of both PPSWs and 
NPSWs for all M c < M < M , then there exist either a PPSW or a PPDL at M = M c 
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if V"'(M c ,0) < 0, whereas for V"'(M c ,0) > 0, there exist either a NPSW or a NPDL at 
M = M c . 

Apart from the conditions of Theorem-6 (Theorem-7), the PPDL (NPDL) solution 
at M — M c is possible only when there exists a PPDL (NPDL) solution in any right 
neighborhood of M c , i.e., PPDL (NPDL) solution at M — M c is possible only when the 
curve M = M D tends to intersect the curve M = M c at some point in the solution space 
of the energy integral, where each point of the curve M = M D corresponds to a PPDL 
(NPDL) solution of the energy integral whenever M D > M c . So, without going through all 
qualitatively different solution spaces of the energy integral by considering the entire range 
of parameters involved in the system, it is not possible to make a systematic investigation 
of the solitary wave solution and/or double layer solution of the energy integral at M — M c . 
The most interesting fact is the existence of the double layer solution at M = M c . 

The existence of solitary structures at M — M c indicates that where and when one can 
apply small amplitude theory to find the weakly nonlinear solitary structures because of the 
following reasons: 

(a) If in a portion of the solution space, there exist solitary waves or double layers of finite 
amplitude at M — M c , then the amplitude of solitary waves corresponding to all M > M c 
is greater than that finite value. Then it is unwise to study weakly nonlinear solitary 
structures by means of any small amplitude methods (viz., replacing V(<f>) by V n (<f>), where 
V n {(j>) represents the power series expansion of the Sagdeev potential or pseudo-potential 
V((f>) in keeping terms up to <p n ) in that portion of the solution space, since these methods 
may lead to erroneous results. Actually, it is not possible to apply small amplitude method to 
study the weakly nonlinear solitary structures in a parameter regime where finite amplitude 
solitary structures exist at M — M c , because in any small amplitude method, amplitude of 
solitary structure always tends to zero as M — > M c + <^> M > M c and M — >■ M c which 
contradicts the occurrence of finite amplitude solitary structure at M — M c . Therefore, 
to study the weakly nonlinear solitary structures by means of small amplitude methods 
for M > M c , the simple prescription is to investigate the nonlinear solitary structures at 
M = M c . So, one cannot apply small amplitude theory to find weakly nonlinear solitary 
structures without investigating the existence of solitary structures at M — M c in the entire 
compositional parameter space. 

(b) If in a portion of the solution space, there exist solitary structures having a finite 
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amplitude at M = M c , then it is simple to find the amplitude as well as width of the 
solitary structure for M = M c + M € (for sufficiently small strictly positive quantity M e ) 
and this can be easily done by plotting V{4>) against <fi at M = M c + M e . However, if in a 
portion of the solution space, there does not exist solitary structures at M — M c , then it is 
very difficult to find the amplitude as well as width of the solitary wave for M = M c + M e 
by plotting V(<f>) against <p. In this situation, the general prescription is to use the small 
amplitude method for finding the weakly nonlinear solitary structure at M — M c + M e . 

(c) Considering any specific problem, it can be easily checked that the solitary structure 
along the curve M = M c is same as the solitary structure for M > M c , if one can move the 
solution space through the family of curves parallel to the curve M = M c . 

Finally, we have used the analytical results as presented in this paper to a specific problem 
of dust acoustic wave in nonthermal plasma to investigate new results and new physical ideas. 

The present paper is organized as follows: Physical interpretation of the energy integral 
for the existence of solitary wave and double layer solutions is given in Sec. [Til In Sec. IHH 
some mathematical preliminaries (including definitions, properties and some basic theorems) 
have been given. Using these definitions, properties and basic theorems we have proved ten 
important theorems in Sec. [TV]for the existence of solitary wave and/or double layer solution 
at M — M c . The analytical results have been verified numerically by considering a problem 
of dust acoustic wave in nonthermal plasma in Sec. |V] Finally, a brief conclusion has been 
given in Sec. IVTJ 

II. PHYSICAL INTERPRETATION OF THE ENERGY INTEGRAL 

The energy integral ([T]) can be regarded as the one-dimensional motion of a particle 
of unit mass whose position is <fi at time £ with velocity d(f)/d£ in a potential well V(0). 
The first term of the energy integral can be regarded as the kinetic energy of a particle 
of unit mass at position <ft and time £ whereas V(4>) is the potential energy of the same 
particle at that instant. Since kinetic energy is always a non- negative quantity, V{4>) < 
for the entire motion, i.e., zero is the maximum value for V(<f>). Again from ([T]), we find 
d 2 (fi/d!; 2 + V'((f>) = 0, i.e. the force acting on the particle at the position <fi is —V'((f>). 
Suppose, V(0) = V'(0) = 0, therefore, the particle is in equilibrium at = because the 
velocity as well as the force acting on the particle at (f> = are simultaneously equal to 
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zero. Now if = can be made an unstable position of equilibrium, the energy integral can 
be interpreted as the motion of an oscillatory particle if V(<f> m ) — for some <p m ^ 0, i.e., 
if the particle is slightly displaced from its unstable position of equilibrium then it moves 
away from its unstable position of equilibrium and it continues its motion until its velocity 
is equal to zero, i.e., until takes the value <f> m . Now the force acting on the particle of 
unit mass at position = m is — V'(4> m ). For <p m < 0, the force acting on the particle at 
the point = m is directed towards the point = if — V(0 m ) > 0, i.e., if V'(0 m ) < 0. 
On the other hand, for m > 0, the force acting on the particle at the point = m 
is directed towards the point = if — V'(0 m ) < 0, i.e., if V'(0 m ) > 0. Therefore, if 
V(4>m) > (for the positive potential side ) or if V'(cj) m ) < (for the negative potential side 
) then the particle reflects back again to = 0. Again, if V(<f> m ) = V'(<p m ) = then the 
velocity d(j)/d£ as well as the force <i 2 0/<i£ 2 both are simultaneously equal to zero at = <f> m . 
Consequently, if the particle is slightly displaced from its unstable position of equilibrium 
(0 = 0) it moves away from = and it continues its motion until the velocity is equal 
to zero, i.e., until takes the value = m . However it cannot be reflected back again at 
= as the velocity and the force acting on the particle at = m vanish simultaneously. 
Actually, if V'(<p m ) > (for m > 0) or if V'(<f> m ) < (for <p m < 0) the particle takes 
an infinite long time to move away from the unstable position of equilibrium. After that 
it continues its motion until takes the value m and again it takes an infinite long time 
to come back its unstable position of equilibrium. Therefore, for the existence of a PPSW 
(NPSW) solution of the energy integral ([I]), we must have the following: (a) = is the 
position of unstable equilibrium of the particle, (b) V(0 m ) = 0, V'(<f> m ) > (V'(0 m ) < 0) 
for some m > (0 m < 0), which is nothing but the condition for oscillation of the particle 
within the interval min{0, m } < < max{0, m } and (c) V(0) < for all < < m 
(0 m < < 0), which is the condition to define the energy integral (PQ) within the interval 
min{0, m } < < max{0, m }. For the existence of a PPDL (NPDL) solution of the energy 
integral (JTJ, the conditions (a) and (c) remain unchanged but here (b) has been modified 
in such a way that the particle cannot be reflected again at = 0, i.e., the condition (b) 
assumes the following form: V((f> m ) = V'(0 m ) = 0, V"(<f) m ) < for some m > (0 m < 0). 

The above discussions for the existence of solitary waves and double layers are valid if 
= is an unstable position of equilibrium, i.e., if V"(0) < along with V(0) = V'(0) = 0. 
In other words, = can be made an unstable position of equilibrium if the potential energy 
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of the said particle attains its maximum value at = 0. Now, the condition V"(0) < gives 
a lower bound M c of M, i.e., V"{0) < <=> M > M c , V"(0) > & M < M c , and 
V"(0) = M = M c <^> V"(M C , 0) = 0. This M c is, in general, a function of the 
parameters involved in the system, or a constant. Therefore, if M < M c , the potential 
energy of the said particle attains its minimum value at = 0, and consequently, = is 
the position of stable equilibrium of the particle, and in this case, it is impossible to make 
any oscillation of the particle even when it is slightly displaced from its position of stable 
equilibrium, and consequently there is no question of existence of solitary waves or double 
layers for M < M c . In other words, for the position of unstable equilibrium of the particle 
at = 0, i.e., for V"(0) < 0(<^ M > M c ), the function V((f>) must be convex within a 
neighborhood of = and in this case both type of solitary waves (negative or positive 
potential) may exist if other conditions are fulfilled. Now suppose that V"(M C , 0) = and 
also V'"(M C , 0) = 0, then if V""(M C , 0) < 0, the potential energy of the said particle attains 
its maximum value at = and consequently, = is the position of unstable equilibrium. 
On the other hand if V"(M c ,0) = 0, V"'{M c ,0) = 0, and V""(M c ,0) > 0, the potential 
energy of the said particle attains its minimum value at = and consequently, = 
is the position of stable equilibrium of the particle and in this case there is no question of 
existence of solitary wave solution and/or double layer solution of the energy integral (0Q). 
But if V"'(M C , 0) ^ along with V(M C , 0) = V'(M C , 0) = V"(M C , 0) = 0, then without 
going through the complete analytical investigation, it is difficult to predict the existence of 
solitary wave and/or double layer solution of the energy integral §Bj at M = M c . In this 
situation, i.e., when V"'{M c ,0) ^ along with V(M C , 0) = V'{M c ,0) = V"(M C , 0) = 0, to 
confirm the existence of solitary wave and/or double layer solution of the energy integral 
(0Q) at M — M c , some definitions and well-known theorems of real analysis are necessary. 
So, in the next section, we consider only those definitions and theorems of real analysis 
for continuation of physical interpretation of the existence of solitary wave and/or double 
layer solution of the energy integral |T| at M = M c when V"'(M C , 0) ^ along with 
V{M C , 0) = V'(M C , 0) = V"(M C , 0) = 0. 
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III. MATHEMATICAL PRELIMINARIES 



A. Definitions and properties 

Dl : Let a be any point of 3ft, the set of all real numbers and 5 > is a strictly positive 
real number, then the S— neighborhood of the point a is denoted by N(a, 5) and is defined 
as a set of all real numbers x such that a — 5 < x < a + 5, i.e., N(a, 5) = (a — 5, a + 5) = 
{x : a — 5 < x < a + 5}, i.e., 

x £ N(a, 5)<=>a-5<x<a + 5<^\x-a\<5. (2) 

PI : If < Si < 5 2 , using Eq. fl2j) it is simple to check the following property: 

N(a,5 1 )CN(a,8 2 ). (3) 

D2 : Let A be any subset of 3ft and a £ A, then a is said to be an interior point of A if there 
exists a 5 > such that N(a, 5) C A, i.e., any point of N(a, S) is also a point of A. 
D3 : 9ft x 3ft is the set of all two-dimensional points (x, y), i.e., 3ft x 3ft is the set of all ordered 
pairs (x,y) such that x £ 3ft and y £ 3ft, i.e., 

3? x 3ft = {(x, y) : x £ 3ft and y £ 3ft}. (4) 

If A and B are any two subsets of 3ft, then we define A x B as the set of all ordered pairs 
(x, y) such that a; £ A and y £ B, i.e., 

A x 5 = {(x, y) : x £ A C 3ft and y £ 5 C 3ft}. (5) 

P2 : If A, B, C and L> are any four subsets of 3ft and A x B C C x D, then it is simple to 
prove the following property: 

ACC and BCD. (6) 

P3 : For any four subsets A, B, C and -D of 3ft, it is simple to prove the following property: 
(A x B) n (C x D) = (A n C) x (5 n D). 

D4 : Let (a, 6) be any point of 3ft x 3ft and 8\ > 0, 5 2 > are two strictly positive real numbers, 
then the SiS 2 - rectangular neighborhood of the point (a, b) is denoted by N(a, b; S±, S 2 ) and is 
defined as a set of all points (x, y) of 3ftx3ft such that a— 8\ < x < a+5± and b—5 2 < y < b+5 2 , 
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i.e., N(a,b;S 1 ,S 2 ) = {(x,y) : x G N(a,Si) and y G N(b,5 2 )} = {(x,y) : a - 5i < x < 
a + 5\ and b — 5 2 < y < b + 5 2 } = {(x, y) : \x — a\ < 5\ and \y — b\ < 5 2 }, i.e., 

(x,y) G N(a,b;5 u S 2 ) 

\x — a\ < Si and \y — b\ < 5 2 . (7) 
P4 : N(a, b; 8 U 5 2 ) = N(a, 8x) x N(b, 5 2 ) 

D5 : Let A x B be any subset of 3? x 3? and (a, 6) be a point of A x B, i.e., a & A and 
b £ B, then (a, 6) is said to be an interior point of A x 5 if there exists a <5i#2- rectangular 
neighborhood of the point (a, b) with 5 X > 0, 5 2 > such that N(a, b; S 1: 5 2 ) C Ax B, i.e., 
any point of N(a,b;S 1: S 2 ) is also a point of A x 5, i.e., N(a,S 1 ) x N(b,5 2 ) C Ax B, i.e., 
7V(aA) C A and JV(6, 5 2 ) C B. 

B. Some theorems of real analysis 

BThl : Let f(x) be a real valued function defined on a subset D of the set of real 
numbers 3? and f(x) be continuous at x = a 6 D. If f(a) ^ 0, then there exists a real 
number 5i > such that f(a)f(x) > for all x G N(a,5i) fl £>. Furthermore, if a is an 
interior point of D, then there exists a real number 5 > such that f(a)f(x) > for all 
x G iV(a,<T) C D. 

BTh2 : Let f(x,y) be a real valued function defined on a subset D = J 4xBofKxSl 
and f(x,y) be continuous at (x, y) = (a,b) G -D. If f(a,b) ^ 0, then there exist two real 
numbers S 1 > and 5 2 > such that f(a,b)f(x,y) > for all (x,y) G N(a : b;S 1 ,S 2 ) fl D. 
Furthermore, if (a, b) is an interior point of D, then there exist two real numbers 63 > and 
S4 > such that /(a, b)f(x, y) > for all (x, y) G iV(a, 6; <5 3 , 5 4 ) C £). 

Bolzano's Theorem : Let f(x) be a real valued function defined on a closed interval [a, b] 
(a < b). If /(x) is continuous on [a, 6] and f(a)f(b) < 0, then there exists a point c G (a, 6) 
such that f(c) =0. 

Lagrange's Mean - Value Theorem : Let f(x) be a real valued function defined on a 
closed interval [a,b] (a < b). If f(x) is continuous on [a, b] and if it is derivable on (a, 6), 
then there exists a point c G (a, 6) such that 



[/(&)-/(«)]/(&- a) = /'(c), 
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where " ' " indicates derivative with respect to x. 

IFT (The implicit function theorem of two variables) : Let V(M, 0) be a real valued 
function of two variables M and defined on Ai x <3>, i.e., for each M G M. and for each 
G $, there exists a unique real number V(M, 0). Let (Mi,0i) G M x $ be an interior 
point of .M x $ and 
Al : \/(M 1 ,0 1 ) = O, 

A2 : |^ is continuous in a rectangular neighborhood of (Mi,0i), 
A3 : f ^ 0, 

there exists a rectangular neighborhood N{M 1 ,(f) 1 \hi,k\) = N(M 1 ,hi) x N(<f>i,ki) of 
(Mi,0i) such that corresponding to every M G N(Mi,h±), 0(g iV(0i,fci)) can be ex- 
pressed uniquely as a function of M, for which the following conditions are simultaneously 
satisfied. 



CI 
C2 
C3 



4>{Mi) = 0i, 

V(M, 0(M)) = for all M G N(Mi, h^, 

l£r = -(m + %) for a11 ( M ' 0) e N ( M ufa h, h) = N{M U h l ) x iV(0 1? 
To prove Theorem 1 - Theorem 10, it is helpful to remember that if a property is true for 
each point of a set then that property is also true for each point of any subset of that set. 
Without any loss of generality, we can assume (M, 0) as an interior point of M x $. 



IV. PROOF OF THEOREM 1 - THEOREM 10 

Theorem 1 : If V"'(M C , 0) ^ 0, then there exists a strictly positive real number e (> 0) 
such that V"'(M c ,0)V"'(M c ,(f)) > for all -0 e < < e . 

Proof : Here /(0) = V'"(M C , 0) is continuous at = and consequently from BThl, we 
can conclude that there exists a strictly positive real number e (> 0) such that 

/(0)/(O)>O V G (O-0 e ,O + e ), 
\/'"(M c ,0)\/'"(M c ,O) >0 V 0G(-0 e ,0 e ), 
=S> V'"(M C , 0)y"'(M c , 0) > V - e < < e . 

(8) 

Theorem 2 : If V"(M C , 0) = and V""(M C , 0) ^ 0, then there exists a strictly posi- 
tive real number e (> 0) such that V"'(M C , 0)V"(M c , 0) > for all < < e and 
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V"'{M C , 0)V"(M c , 0) < for all -0 e < < 0. 

Proof : From Theorem 1, we see that if V"'(M Cl 0) ^ 0, then there exists a strictly positive 
real number e (> 0) such that 

V'"(M e , <f>)V'"(M c , 0) > V - e < < (j) e . (9) 

Again as V(M C , 0) and its derivatives of all order with respect to exist finitely, /(0) = 

V"(M C , 0) is a continuously derivable function of G $. Now, for any 0(^ 0) G (— e , e ), 

we have (a) the real valued function /(0) is well-defined in the closed interval [min{O,0}, 

max{0, 0}], (b) the real valued function /(</>) is continuous in the closed interval [min{0, 0}, 

max{0, 0}], (c) the real valued function /(</>) is derivable in the open interval (min{O,0}, 

max{O,0}). Therefore, from the Lagrange's mean - value theorem of differential calculus, 

we can conclude that there exists a 0i G (min{0, <f)}, max{0, 0}) such that 
/(max{0, (j)}) - /(min{0, 0}) = ^ 
max{0, 0} — min{0, 0} 

^Mzm = mzm =fM . (10) 

As /(0) = V"(M C , 0) = and ^ 0, (HUD assumes the following form: 



=> y"(M c ,0) = 0\/" / (M c ,0 1 ), 

^ V'"(M C , 0)V"(M c , 0) = <f>V"'{M c , 0)V"'(M c , 00. (11) 

Suppose that < 0, then min{0, 0} = and max{0, 0} = 0, and consequently, min{0, 0} < 
0x < max{0, 0} =>■ < 0i < 0, but since 0(^ 0) G (— e , e ) we have — e < < e =>• 

-0 e < < 01 < < £ -0 e < 01 < e 01 G (-0 e , e ) . 

Next suppose that > 0, then min{0, 0} = and max{0, 0} = 0, and consequently, 
min{O,0} < 0i < max{O,0} =>- < 0i < 0, but since 0(^ 0) G (— e ,0 e ) we have 

-0 e < < £ ^ -0 e < < 01 < < e ^ -0 e < 01 < e ^ 01 G (-0 e , e ). 

Therefore in any case (0 < or > 0), we always have 0i G (— e , e ), and consequently, 
from Eq. (jUJ), we get 

r /,/ (M c ,0i)V /,/ (M c ,O) >0. (12) 
From Eq. (JT2"|) . we get 

01/ W (M C , 0i)V"'(M c , 0) < 0V - e < < 0, (13) 
0\/'"(M c , 0i)V r/// (M c , 0) > 0V < < e . (14) 
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From Eq. (TTTT) . using Eqs. (fl3|) and (TH|) . we get 

niJ)HMJ)<OV -0 e <0<O, (15) 
y'"(M c , 0)V"(M c , 0) > V < < e , (16) 

and consequently, Theorem 2 is proved. 

From Theorem 2, we observe that whenever V(M, 0) = V'(M, 0) = ^"(Mc, 0) = and 
V"'(M c ,0) < 0, then V"(M e ,<j>) > for all -0 e < < and K"(M c ,0) < for all 
< < <f) e , i.e., V(M C , 0) is a concave function within the interval — e < < and 
V(M C , 0) is a convex function within the interval < < e , where e is strictly positive 
quantity but it can be made sufficiently small. Therefore, if the particle placed at = be 
slightly displaced towards the positive potential side, it falls within the interval < < e 
and due to the convexity of the function V(M C , 0) within the interval < < e , it moves 
away from = and it continues its motion until its velocity is equal to zero, i.e., until 
takes the value = C > 0, where V(M c ,(f) c ) = and in this case if V'(M c ,<j) c ) > 
(V'{M c ,(j) c ) = and V"(M C , <p c ) < 0), one can easily get a PPSW (PPDL) as a solution 
of the energy integral ([T|). Again, as V(M C , 0) is a concave function within the interval 
— e < < 0, if the particle placed at = be slightly displaced towards the negative 
potential side, it falls within the interval — e < < and due to the concavity of the 
function V(M C , 0) within the interval — e < < 0, it moves towards the point = and 
consequently, there does not exist any solitary wave or double layer solution in the negative 
potential side. In true physical sense, lim^o+o V(M C , 0) is an unstable position which lies 
in a small right neighborhood of the equilibrium position = 0, whereas lim^o-o V{M C , 0) 
is a stable position which lies in a small left neighborhood of the equilibrium position = 0, 
where 0— i-O + O<^0>O and approaches to zero, i.e., approaches to zero from the right 
side of zero and 0— t-0 — O<^0<O and approaches to zero, i.e., approaches to zero from 
the left side of zero. Hence, as lim^o+o V(M C , 0) is an unstable position which lies in a small 
right neighborhood of the equilibrium position = 0, there may exist either a solitary wave 
or a double layer in the positive potential side. On the other hand, as lim^^ _ V(M C , 0) is 
a stable position which lies in a small left neighborhood of the equilibrium position = 0, 
there does not exist any solitary wave or double layer in the negative potential side. Hence 
the above theoretical discussions can be summarized as follows. 

Theorem 3 : If V(M, 0) = V'(M, 0) = V"(M C , 0) = and V"'(M C , 0) < 0, then at M = M c , 
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there may exist either a solitary wave or a double layer in the positive potential side but 
there does not exist any solitary wave or double layer in the negative potential side. 

Again from Theorem 2, we observe that whenever V(M, 0) = V'(M, 0) = V"(M C , 0) = 
and V"'(M c ,0) > 0, then V"(M c ,(j)) < for all -<f> e < < and V"(M c ,</>) > for all 
< <fi < <f) e , i.e., V(M C , <f>) is a convex function within the interval — <f) e < <ft < and V(M C , <ft) 
is a concave function within the interval < < <p t . Therefore, following the same argument 
as given in the proof of Theorem 3, we have 

Theorem 4 : If V(M, 0) = V'{M, 0) = V"(M C , 0) = and V"'(M C , 0) > 0, then at M = M c , 
there may exist either a solitary wave or a double layer in the negative potential side but 
there does not exist any solitary wave or double layer in the positive potential side. 

From Theorem 3 and Theorem 4, we note that if V(0) = V'(0) = 0, V"(M C , 0) = and 
V"'(M c ,0) 7^ 0, the point <fi = is the point of inflexion that seperates the convex part of 
V(M C , <fi) from the concave part and if <fi = is a point of inflexion then there may exist 
either a solitary wave or a double layer in either positive potential side or negative potential 
side, but the solitary structures (including double layers) of both polarities can not coexist 
at M = M c . Therefore, 

Theorem 5 : If V(M, 0) = V'(M, 0) = V"(M C , 0) = and V'"(M C , 0) ^ 0, it is not possible 
to have coexistence of both positive and negative potential solitary structures (including 
double layers) at M — M c . 

Next to confirm the existence of the solitary wave and/or double layer solution of the 
energy integral ([T|) at M = M c , we have proved the following theorems. 
Theorem 6 : If V(M, 0) = V'(M, 0) = V"(M C , 0) = 0, V"'(M c ,0) < 0, dV/dM < for 
all M > and for all <fi > 0, and if there exists at least one value M of M such that the 
system supports Positive Potential Solitary Waves (PPSWs) for all M c < M < M , then 
there exist either a PPSW or a Positive Potential Double Layer (PPDL) at M = M c . 
Proof : From Theorem 3, we have seen that if V(M, 0) = V'(M, 0) = V" (M c , 0) = 
and V"'(M C , 0) < 0, there may exist either a solitary wave or a double layer in the positive 
potential side at M = M c but there does not exist any solitary wave or double layer in 
the negative potential side at M = M c . Now with the help of IFT, we shall confirm the 
existence of either a solitary wave or a double layer in the positive potential side at M = M c . 
Let M c < Mi < Mo, then according to the condition of the present theorem, there exists a 
solitary wave in the positive potential side at M = M x . Let the amplitude of this solitary 
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>0, 



wave is <f>i(> 0), then from the conditions for existence of PPSW we have 
SM1-1 :: V(M u fa) = 0, 
SM1-2 :: g 

SM1-3 :: V(M 1 , fa < for all < <f) < fa. 

Again, as dV/d(j) is a continuous function of (M, fa, we get 

SM1-4 :: ^ is continuous on a rectangular neighborhood of (Mi, fa). 

Therefore, SM1-1, SM1-4, and SM1-2 are, respectively, same as the conditions Al, A2, 
and A3 of IFT, and consequently, we can apply IFT in a rectangular neighborhood of 
the point (Mi, fa). Therefore, from IFT, we can conclude that there exists a rectangular 
neighborhood N(Mi, fa; hi, ki) = A r (M 1 , hi) x N(fa, ki) of (Mi, 4>i) such that corresponding 
to every M G N(M i: hi), fae N(fa, ki)) can be expressed uniquely as a function of M, for 
which the following conditions are simultaneously satisfied. 

IFT-C1 :: faMi) = fa, 

IFT-C2 :: V(M, faM)) = for all M G N(M u hi), 

IFT-C3 :: $ = -(f£ - %) for all (M, fa G N(M U fa; hi, h) = N(M U hi) x N(fa, h). 

Again from SM1-2, using theorem BTh2, we have strictly positive real numbers /i2(> 0) 
and k 2 (> 0) such that 

SM1-5 :: g > for all (M, fa G N(M U fa; h 2 , k 2 ) = N(M U h 2 ) x N(fa, k 2 ). 

Again from IFT-C3, we see that 4> is derivable function of M for M G (Mi— hi, Mi+hi), and 
consequently, it is continuous in (M 1 — hi,Mi + hi). As is continuous in (Mi —hi,Mi + hi) 
and faMi) — fa > 0, using theorem BThl, we have a strictly positive real number h 3 (> 0) 
such that 

SM1-6 :: faM) > for all M G N(M u h 3 ). 

So, if we take h = min{/ix, h 2 , h 3 } and k = minj/cx, k 2 }, then N(Mi, h) C N(M 1 , hi) for 
all i = 1,2,3, and N(fa,k) C N(fa,kj) for all j = 1,2, consequently, N(Mi, fa; h, k) = 
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N(M 1 , h) x N((j> u k) C N(M u hi) x N^kj) = N{M u <p 1 ;h u k j ) for all i = 1,2,3 and 
j = 1,2. Therefore, we can also apply IFT in the rectangular neighborhood N(M\, 0i; h, k) = 
N(M 1 ,h) x N((j)i,k), and consequently, with respect to the rectangular neighborhood 
N(M 1 ,(p 1 ;h,k), IFT-C1, IFT-C2, IFT-C3, SM1-5, and SM1-6 can be put in the fol- 
lowing form: 

mSMl-1 :: 0(Mi) = 0i, 

mSMl-2 :: V(M, 0(M)) = for all M G N(M U h), 

mSMl-3 :: $ = ~(m - %) for a11 ( M ' ^) e W(M ls K k) = N(M U h) x N(<j> u k), 

mSMl-4 :: f£ > for all (M, 0) G iV(Mi, i; /i, fc) = iV(Mi, /i) x iV(0i, fc), 

mSMl-5 :: 0(M) > for all M G N(M 1: h). 

Now choose two strictly positive real numbers e\ and e such that M\ = M c + e± and 
N(M c ,e) C N(M 1 ,h). These two conditions are simultaneously satisfied if we choose < 
ei < /i and < e < h - e 1 . As 0(M) G iV(0i, k) for every M G N(M U h) and as M c G 
iV(Mi, /i), then we must have 0(M C ) G AT(0 1; k) <fi c G iV(0i, fc), where we set C = 4>{M C ). 
Now choose two strictly positive real number S\ and 5 such that 0i = c + <5i and N(<f) c , S) C 
N((f>i,k). These two conditions are simultaneously satisfied if we choose < 5\ < k and 
< 5 < k-5x. 

Therefore, N(M c ,e) C N(M 1 , h) and iV(0 c , 5) C N(fa,k) => N(M c ,e) x iV(0 c , 5) C 
N(M 1: h) x N((j)i, k) => N(M C) 4> c ] e, <5) C N(M 1: X ; /i, fc), and consequently, we can also ap- 
ply IFT in the rectangular neighborhood N(M C , C ; e, 5) C N(M 1: X ; /i, fc). So, with respect 
to the rectangular neighborhood iV(M c , C ; e, 5), IFT-C1, IFT-C2, IFT-C3, mSMl-4, 
and mSMl-5 can be put in the following form: 

fmSMl-1 :: 0(M C ) = C , 

fmSMl-2 :: V{M,<j>{M)) = for all M G JV(M c ,e), 

fmSMl-3 :: $ = ~(m - f ) for a11 ( M - ^) e ^( M c, <fe e, 5) = iV(M c , e) x iV(0 c , 5), 
fmSMl-4 :: g > for all (M, 0) G iV(M c , C ; e, 5) = iV(M c , e) x iV(0 c , 5), 
fmSMl-5 :: <j>{M) > for all M G JV(M c ,e). 

15 



As M c G N(M C , e) and C G iV(0 c , 5), from fmSMl-1, fmSMl-2, fmSMl-4, and fmSMl- 
5, we get the following conditions: 



SMc-1 :: V(M C 



SMc-2 :: 



0, where C = 0(M C 



> 0, 



(M C ,0 C ) 

SMc-3 :: C = <f>(M c ) > 0. 

It is important to note that if we take M G (M c — e,M c ), then also the above three 
conditions are simultaneously satisfied but M < M c <^ V"(M, 0) > 0, and consequently, 
= is a position of stable equilibrium. So, there is no question of existence of any solitary 
wave or double layer solution of the energy integral ([I]) for all M G (M c — e, M c ). 

Again, from the condition of the present theorem, we have dV/dM < for all M > 
and for all > 0. As our discussion is restricted to positive potential side (0 > 0) only, 
from fmSMl-3 and fmSMl-4, we get 

Mc :: ^ > for all M G N(M C , e), 

which shows that is an strictly increasing function of M for all M c — e < M < M c + e. 
Now we shall prove that 

SMc-4 :: V(M C , 0) < for all < < C , 

and the following analysis have been developed to prove SMc-4. 

Consider a value M 2 of M such that M c < M 2 < M c + e, then there exists a PPSW of 
amplitude 02, and consequently we have 

SM2-1 :: V(M 2}( j) 2 ) = 0, 



SM2-2 :: 



dv 



> 0, 



(M 2 ,0 2 ) 

SM2-3 :: V{M 2 , 0) < for all < < 2 . 

As M c < M 2 < M c + e =^ M c - e < M 2 < M c + e, we have 2 = 0(M 2 ). Again as is 
an increasing function of M for all M c - e < M < M c + e, and M c , M 2 G N(M C , e) with 
M c < M 2 , we have 0(M C ) < 0(M 2 ) =4> C < 2 . we can write SM2-3 as follows 



mSM2-3 :: V(M 2 , 0) < for all < < C < 



>2- 
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If possible, let there exists a value P of such that < <p p < C and 
Ql :: V(M c ,(f> p ) > 0. 

Again from mSM2-3, as < P < C < 02, we get 
Q2 :: V(M 2 ,(j) p ) < 0. 

Considering V(M, P ) as a function of M only, from Ql, Q2 and Bolzano's Theorem, we 
can conclude that there exists a value M 3 of M such that 

V(M 3 ,<f> p ) = 0, (17) 

where M c < M 3 < M 2 < M c + e. But as M c < M 3 < M 2 < M c + e, there exists a PPSW at 
M = M 3 of amplitude 3 > 0, and consequently we have 

SM3-1 :: V(M 3 ,<j> 3 ) = 0, 



SM3-2 :: ^ 



> 



(M 3 ,<p 3 ) 

SM3-3 :: V{M 3 , 0) < for all < < 3 , 

As M c < M 3 < M c + e =^ M c - e < M 3 < M c + e, we have 3 = 0(M 3 ). Again as is an 
increasing function of M for M c - e < M < M c + e, and M c , M 3 G A^(M C , e) with M c < M 3 , 
we have 0(M C ) < 0(M 3 ) =^ C < 3 . So we can write SM3-3 as follows 

mSM3-3 :: V{M 3 , 0) < for all < < C < 3 , 

Again from mSM3-3, as < P < C < 3 , we get 

Q3 :: V(M 3 , p ) < 0, 

which contradicts Eq. (TTTl) . i.e., V(M 3 ,<f> p ) = 0. Therefore, our supposition is wrong, 
and consequently, there does not exist any P 6 (O,0 C ) such that V(M c ,(p p ) > 0, i.e., 
l^(M c , P ) ^ and consequently, we get 

mSMc-4 :: V(M C , 0) < for all < < C . 

Now mSMc-4 suggests the following two possible cases: 

Case - 1: There does not exist any e (O,0 C ) such that V{M Cl 0) = 0, i.e., V^Mc^) < 
for all < < C . 
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Case - 2: There exists at least one G (0, C ) such that V(M C , 0) = 0. 
Combining Case - 1 with SMc-1, SMc-2 and SMc-3, we get 

(1) : V(M c ,<j> e ) = 0, where C = 0(M C ), 

(2) = % >0, 

(M c ,<f> c ) 

(3) : C = 0(M C ) >0, 

(4) : V{M C , 0) < for all < < C . 

These four conditions confirm the existence of PPSW solution of the energy integral (pQ) at 
M = M c . 

For Case - 2, there exists at least one G (0, C ) such that V(M C , 0) = 0. Let pl be 
the smallest strictly positive root of the equation V(M c ,(p) = such that p i G (O,0 C )(<^ 
< p i < C ) and consequently, for this case, we have the following subcases. 

Subcase-1: V(M C ,4>) < for all < < pi and V(M C , 0) = for all G [0 p i,0 c ](See 
FIG.^a)). 

Subcase-2: V(M c ,<f)) < for all < < pi and V(M C , 0) = for all G [0 p i,0 P 2] with 
< c (See FIG. 1Kb)). 



Subcase-3: V(M C , 0) < for all < < pi and F(M C , 0) < for all G (0 pi - e pi , <p pl ) 
as well as G (0 p i,0 p i + e p i) along with l^(M c ,0 pl ) = for some e p i > (See FIG. 

mc)). 

As V(M c ,(p) is a continuously differentiable function of for all G $ and for the 
Subcase-1, it is simple to check that V(M C , 0) is not differentiable at = p i, Subcase-1 
is not possible. With the help of same argument, it is simple to check that Subcase-2 is 
not possible when pl < p2 . When pl = p2 , then either we have Subcase-3 or we have 
two different branches of the curve V(M C , 0) originated from = p i. The later case is 
not possible because V(M c ,<fr) is not differentiable at = pl . So, Subcase-3 is the only 
possible case of Case-2 and for this subcase we have the following conclusions. 

Subcase-3 : In this V(M C , 0) is a continuously differentiable function of 0, 

V(M c ,(f)) is not only attains its maximum value at = pl , V(M C , 0) is convex for all 
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and w 



G (0 p i — e p i, (fipi + e p \). Therefore, from the elementary theory of convexity of a function, 
we have 
dV_ 

^ (Af ,^pi) 

Combining these two conditions along with the first condition of Subcase-3, i.e., V(M C , 0) < 
for all < < (f> p i, we get 

(1): V(M c ,4> pl )=0, 



< 0. (18) 

(M c ,if>pi) 



(2): S 



= 0, 



<0, 

(M c ,(j> p i) 



( 3 ) : 9 d( p 

(4) : pl >O, 

(5) : V{M C , 0) < for all < < <j> pX . 

These five conditions confirm the existence of PPDL solution of the energy integral ([1]) at 
M = M c . Therefore, under the conditions of the present theorem, one can get either a 
PPSW or a PPDL as a solution of the energy integral (pQ) at M — M c . 
Theorem 7 : If V(M, 0) = V'(M, 0) = ^"(M c , 0) = 0, V"'(M c ,0) > 0, S^/^M < for 
all M > and for all <fi < 0, and if there exists at least one value M of M such that the 
system supports Negative Potential Solitary Waves (NPSWs) for all M c < M < M , then 
there exist either a NPSW or a Negative Potential Double Layer (NPDL) at M = M c . 
Proof : Same as Theorem 6 with a slight modification for the case of negative potential 
side. 

Theorem 8 : If V(M, 0) = V'(M, 0) = V"{M C , 0) = 0, V"'(M c ,0) < 0, dV/dM < for all 
M > and for all > 0, and if the system supports only NPSWs for M > M c , then there 
does not exist any PPSW at M = M c . 

Proof : From Theorem 3, we have seen that if V(M, 0) = V'(M, 0) = V"(M C , 0) = and 
V"'(M C , 0) < 0, there may exist solitary wave in the positive potential side at M = M c but 
there does not exist any solitary wave in the negative potential side at M = M c . If possible 
let there exists a PPSW at M — M c and for M > M c , the system supports only NPSWs. 
Therefore, if C is the amplitude of the PPSW at M = M c , the following conditions are 
satisfied: 



19 



Mcl:: V(M c ,0 c ) = 0, 
Mc2:: C > 0, 



Mc3:: ^ 



>0, 

(M c ,<t> c ) 



Mc4:: V(M C , 0) < for all < < C . 

As dV/dcj) is a continuous function in its domain of definition, we have, 

Mc5:: |^ is continuous in some rectangular neighborhood of (M C ,0 C ). 

Therefore, all the assumptions, viz., Mcl, Mc5 and Mc3, of IFT at (M c , C ) are sat- 
isfied by the function V(M,<f>). Therefore, from IFT, we can conclude that there exists a 
rectangular neighborhood N(M C , C ; Si, 62) = N(M C , Si) x iV(0 c , 82) of (M c , C ) such that cor- 
responding to every M G N(M C) 81), 0(g iV(0 c , <5 2 )) can be expressed uniquely as a function 
of M, for which the following conditions are simultaneously satisfied. 

IFT Mcl :: 0(M C ) = C , 

IFT Mc2 :: V(M,<f>(M)) = for all M G N(M c ,<?i), 

IFT Mc3 :: ^ = -(Jg - §£) for all (M, 0) G iV(M c , C ; ^ <J 2 ) = N(M C , Si) x iV(0 c , <f 2 ). 

Again from Mc3, using BTh2, we have two strictly positive real numbers 8 3 (> 0) and 
&t(> 0) such that 

Mc6 :: % > for all (M, 0) G iV(M c , C ; <5 3 , <5 4 ) = iV(M c , S 3 ) x iV(0 c , <S 4 ). 

So, if we take 5 = mm{S 1 ,8 3 } and p = min{5 2 ,<5 4 }, then N(M C) S) C N(M c ,8i) for all 
? = 1,3, and iV(0 c , p) C iV(0 c , ^) for all j = 2, 4, consequently, N(M C , C ; 5, p) = N(M C , S) x 
A^(0 C , p) C N(M C , Si) x iV(0 c , Sj) = N(M C , C ; Sj) for alH = 1, 3 and j = 2, 4. Therefore, 
we can also apply IFT in the rectangular neighborhood N(M C , C ; 8, p) = N{M C , 8)xN((f) c , p), 
and consequently, IFT Mcl, IFT Mc2, IFT Mc3, and Mc6 can be re-written as follows: 

mMcl :: 0(M C ) = C , 

mMc2 :: V{M,<j>(M)) = for all M G N(M C ,6), 

mMc3 :: ^ = -(J£ - f ) for all (M, 0) G iV(M c , C ; 5, p) = 7V(M C , 5) x iV(0 c ,p), 
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mMc6 :: > for all (M, 0) G iV(M c , C ; 5, p) = iV(M c , 5) x iV(0 c , p). 

Again, from the condition of the present theorem, we have OV/dM < for all M > 
and for all > 0. As our discussion is restricted to positive potential side (0 > 0) only, from 
mMc3 and mMc6, we get <t>(M) is an increasing function of M for all M c — 5 < M < M c + 5. 
Consider a value Mi of M such that M c < M x < M c + 5 => M c - 5 < M c < Mi < M c + 5 
M c - 5 < Mi < M c + 5 =^ Mi G N(M C ,8), and consequently, from mMc2, the condition 
that the function 0(M) is an increasing function of M for all M c — 5 < M < M c + 5 and 
mMc6, we get the following three conditions: 

Ml-1:: V(M X M = 0, with <f> x = <j>(Mx), 

Ml-2:: M c - 5 < M c < Mi < M c + 5 M c , Mi G A^(M C , 5) with M c < Mi => 0(M C ) < 
0(Mi) =^ C < 0i 0i > C > 0, 



[It is important to note that if we take M G (M c — 5, M c ), then also the above three conditions 
are simultaneously satisfied but M < M c V"(M, 0) > 0, and consequently, = is a 
position of stable equilibrium. So, there is no question of existence of any solitary wave 
solution of the energy integral (CQ) for all M G (M c — 8, M c ).] 

Now if V(M 1 ,<f>) < for all < < 0i, then from Ml-1, Ml-2, and Ml-3, we can 
conclude that there exists a PPSW at M = M x > M c but this is impossible as the system 
supports only NPSWs for M > M c . Therefore, there exists at least one value 02 of such 
that 



Now we shall prove that actually 02 is restricted by the following inequality: C < 02 < 0i- 
Now from the condition dV/dM < for all M > and for all > 0, we have V(M, 0) 
is a strictly decreasing function of M for all M > and for any given value of > 0. As 

Mi > M c , we get 

Ml-4 :: V(M h </>) < V(M e ,<f>) for all > 0. 

Again, from Mc4, we get V(M C , 0) < for all < < C => V(M U 0) < V(M C , 0) < for 
all < < C , where we have used the inequality (V(Mi,<f>) < V(M c ,(f)) for all > 0) as 
given in Ml-4. Therefore, we have 



Ml-3:: % 



>0, 



(M lt 4>i) 



V(M U 2 ) > for < 02 < 0i. 



(19) 
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Ml-5 :: V(M X , 0) < for all < < <f> c . 

Again, as C > and V(M C , (f> c ) = 0, putting = C in Ml-4, we get 
Ml-6 :: V(M X , C ) < V(M c ,<j> e ) = V{M x ,cj) c ) < 0. 
Combining Ml-5 and Ml-6, we get 
Ml-7 :: V{M X , 0) < for all < < C . 

Therefore, from Ml-7, we can conclude that C < 02 < 0i. Now we shall prove that 
there exists a value M 2 of M such that 2 = 0(M 2 ) for M c < M 2 < M x . 

Let us define the function ip(M) = 0(M) - 2 for all M G [M C ,M X ). Now as 0(M) is 
a continuously differentiable function of M for all M G N(M C ,S), then it is automatically 
continuous for all M G N(M C ,5). A gain as [M c , M x ] C N[M C ,5), 4>{M) is also continuous 
for all M G [M c ,Mi], and consequently, ip(M) is continuous for all M G [M c ,Mi\. Now as 
C = 0(M C ) and 0i = 0(Mi), the condition C < 2 < 0i gives 0(M C ) < 2 < 0(Mi) <^> 
0(M C ) - 2 < and 0(M X ) - 2 > <=> ip{M c ) < and ij){M x ) > 0. 

Therefore, we have a continuous function ijj(M) : [M c , Mi] i — > 3? such that ip{M c ) < 
and ip(Mi) > 0,and consequently, by Bolzano's theorem, we can conclude that there exists 
a value M 2 of M such that ip(M 2 ) = 0o 0(M 2 ) = 2 for M c < M 2 <M X . Now as V(M, 0) 
is strictly decreasing function of M > for any fixed value of > 0, and consequently, as 
M 2 < M x , we have V(M 2 , 0) > V(M X , 0) for any given value of > 0. If we take = 2 , 
we get K(M 2 ,0 2 ) > V{M x ,(f> 2 ) V{M Xl (j) 2 ) < ^(M 2 ,0 2 ) = F(M 2 ,0(M 2 )) = because 
M c < M 2 < M x < M c + 5 => M 2 G N(M C ,5) => V(M 2 , <f)(M 2 )) = 0, where we have used 
the condition mMc2 to get V(M 2 ,(j)(M 2 )) = 0, So, we have V(M x ,(f) 2 ) < which again 
contradicts Eq. (ITUl) . i.e., V r (M 1 ,0 2 ) > 0. Thus, our original assumption is wrong, and 
consequently, theorem is proved, i.e., if the system supports only NPSWs for M > M c , then 
there does not exist any PPSW at M = M c provided the conditions of the present theorem 
are fulfilled. 

Theorem 9 : If V(M, 0) = V'(M,0) = V"(M c ,0) = 0, V"'(M c ,0) > 0, dV/dM < for all 
M > and for all < 0, and if the system supports only PPSWs for M > M c , then there 
does not exist any NPSW at M = M c . 

Proof : Same as Theorem 8 with a slight modification for the negative potential side. 
Theorem 10 : If V(M, 0) = V'(M, 0) = V"(M C , 0) = 0, V"'(M C , 0) ^ 0, dV/dM < for all 
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M > and for all <fi 7^ 0, and if there exists at least one value Mo of M such that the system 
supports both positive and NPSWs for all M c < M < M , i.e., if the system supports the 
coexistence of both PPSWs and NPSWs for all M c < M < M , then if V"'(M C , 0) < 0, there 
exist either a PPSW or a PPDL at M = M c whereas if V"'(M C , 0) > there exist either a 
NPSW or a NPDL at M = M c . 

Proof : Follows from Theorem 6 and Theorem 7. 

From theorem 6 -theorem 10, it is simple to see that the critical Mach number alone 
cannot determine the existence of particular nonlinear structure and its polarity. This is 
not even the motivation of our present manuscript. However, our intension is to determine 
the conditions under which nonlinear structure exists at the critical Mach number. The 
existence of nonlinear structure at the critical Mach number and its polarity is determined 
by the nature of the solitary structures for at least one right neighbourhood of the point 
M = M c together with the sign of the function V"'(M C , 0). 

Except Theorem 1 and Theorem 2, the proof of rest of the theorems (Theorem 3 - Theorem 
10) carry physical sense inherently, since the proofs are based on the physical interpretation 
of the energy integral as discussed in Sec. II. For example, Theorem 6 is proved with the help 
of Theorem 3 while Theorem 3 is proved on the basis of physical interpretation of the energy 
integral. One can get new result behind the nonlinear structures at M — M c by considering 
any specific problem, only when these theorems have been used on the qualitatively different 
solution spaces or compositional parameter spaces with respect to any parameter of the 
system (keeping arbitrary fixed values of other parameters in their respective physically 
admissible range) showing the existence region of different types of solitary structures for 
M > M c . For this purpose, we have considered the same problem of dust acoustic wave in 
nonthermal plasma of Das et al— in Sec. |V] 

So, without going through all qualitatively different solution spaces of the energy integral 
by considering the entire range of parameters involved in the system, it is not possible 
to make a systematic investigation of the solitary structures at M — M c . To construct the 
solution space or compositional parameter space with respect to any parameter of the system 
(keeping arbitrary fixed values of other parameters in their respective physically admissible 
range) showing the existence region of different types of solitary structures for M > M c , it 
is not necessary to know the existence of solitary structure at M = M c . For M > M c , using 
the physical interpretation as given in first para of Sec. II, one can easily set up a numerical 
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scheme to construct the compositional parameter space with respect to any parameter of 
the system (keeping arbitrary fixed values of other parameters in their respective physically 
admissible range) showing the existence region of different types of solitary structures. So, 
to construct the solution space, the existence of solitary structure at M — M c does not play 
any role. 



V. APPLICATION 

Here we consider the specific problem of Das et alr§- to apply the analytical results (theo- 
rem 3 - theorem 10) as presented in this paper. In this paper,— a computational scheme has 
been developed to study the arbitrary amplitude dust acoustic solitary waves and double 
layers in a nonthermal plasma consisting of negatively charged dust grains, nonthermal ions 
and isothermal electrons including the effect of dust temperature. The Sagdeev potential 
approach, which is valid to study the arbitrary amplitude solitary waves and double layers, 
has been employed. Four basic parameters of the system are /i, a, (3% and ad, which are re- 
spectively the ratio of unperturbed number density of electrons to that of nonthermal ions, 
the ratio of average temperature of nonthermal ions to that of isothermal electrons, a pa- 
rameter of the nonthermal distribution of ions and the ratio of average temperature of dust 
particles to that of ions divided by the number of negative charges residing on a dust grain 
surface. The Sagdeev potential V((p)(= V(M, 4>)) as given by Eq. (21) of Das et ai.,— is well 
defined as a real number for all M G Ai and for all G where Ai is the set of all stictly 
positive real number and $ = {(f) : ^> m < <ft < +°°}) an d ^ m is given by the first equation of 
(17) of Das et ai.— So, except the point <fr — ^m, any G $ is an interior point of $ and any 
M > is an interior point of Ai. The solution spaces of the energy integral (20) with respect 
to Pi has been shown in FIG. 3 of Das et ai.— for fixed values of other parameters. From 
this solution space, they have found that for any fixed values of the parameters /z, a and 
the entire interval of Pi can be broken up into four disjoint subintervals I : < fi\ < Pi c , II 
: Pic <Pi< 02c, HI : P 2c <Pi< P 3c , IV : /3 3c < Pi < Pm, where P M = min{l + apt, 4/3}. 
FIG. 3 of Das et ai.— also shows that (i) in subinterval I, only NPSWs can exist and the 
Mach number M for these waves lies within the interval M c < M < M max , where M c is 
the lower bound of M and M max is the upper bound of M which is defined only when the 
system can support NPSWs, (ii) in subinterval II, both NPSWs and PPSWs can coexist and 

24 



the Mach number M for these waves lies within the interval M c < M < Md, whereas only 
NPSWs are possible if M D < M < M max , (iii) in subinterval III, both NPSWs and PPSWs 
can coexist and Mach number for these waves lies within M c < M < M max , whereas only 
PPSWs are possible if M max < M < Md, (iv) in subinterval IV, only PPSWs can exist and 
Mach number M for these waves lies within M c < M < Md- In subintervals II, III and IV, 
only PPDL solution is possible along the curve M = Md- Again from this figure, we see that 
the curve M = Md tends to intersect the curve M = M c at 0\ = 0\ c , i.e., there always exists 
a PPDL solution in any right neighborhood of M c for some 0\ > 0\ c . Therefore, at 0\ = 0i c , 
we can expect a PPDL solution at M — M c , if the other conditions of Theorem 6 have been 
fulfilled. Before going to discuss the existence of solitary wave and/or double layer solutions 
of the present problem at M = M c , it is necessary to know the sign of dV/dM, V'"(M c ,0) 
and the sign of V""(M C , 0) when V"'(M C , 0) = 0. From the expression of V(4>)(= V(M, 0)) 
as given by Eq. (21) of Das et al.,— we get the following equations. 

0] -M(V^-^LY, (20) 



dM I V™d 



V'"(M C , 0) = _12^3+3^-/4(l-«V) > (21) 

where 

/3 T = l + o:/x-/3i,/x T = l-/x. (22) 

From Eq. (120|) . we see that dV/dM < for all M > and for all ^ 0, and consequently, 
for all M > 0, the condition dV/dM < is satisfied for all <ft > as well as <p < 0. 

Now from the expression of V"'(M C , 0) as given by Eq. ( 12 ip of present paper, it can be 
easily checked that (d/d(3i)(V"'(M c , 0)) > for all < /?i < /3 M and for any admissible 
values of the other parameters, consequently, V"'(M C , 0) is strictly increasing function of fii 
for all < j3\ < f3 M . Again, it can be easily checked that V"'(M C , 0) < for 0i = and 
V"'(M C , 0) > for 0i = 1 + afi > (3m- Therefore, V"'(M C , 0) strictly increases with increasing 
0i starting from a negative value and ending with a positive value. So, V'"(M C , 0) intersects 
the axis of 0\ at 0\ = /3 C (0 < (3 C < 0m) and consequently, we have V"'(M c ,0) < for all 
< 0i < 0c, V"'{M C , 0) > for all C < X < M and V"'(M C , 0) = at 0x = C . It is also 
interesting to note that the equation V'"(M C , 0) = has one and only one real root 0\ = C 
of 0i within the admissible range of 0\ and for any fixed values of the other parameters. 
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Again, it is a simple task to check that V""(M C , 0) > at X = C , i.e., V""{M C , 0) > when 
V'"(M C , 0) = 0. Therefore, when V"'(M C , 0) = 0, <ft — is the position of stable equilibrium 
and consequently, there does not exist any solitary structure of the energy integral (20) of 
Das et at M = M c . Now we are in a position to discuss the existence of solitary wave 
and/or double layer solution of the energy integral (20) of Das et al,— at M = M c for any 
admissible value of 0i and fixed values of the other parameters. 

In subinterval I (0 < 0\ < 0i c ), one can expect only NPSW at M = M c , provided 
V"'(M C , 0) > (from Theorem 4). But it can be easily checked that V'"(M C , 0) < for any 
value of 0i lying within subinterval I. So from Theorem 8, we can conclude that there does 
not exist any solitary wave solution at M = M c for any value of 0i lying within subinterval I. 
Again, in subinterval IV (0 3c < 0i < (3m), one can expect only PPSW at M = M c , provided 
V"'(M C ,0) < (from Theorem 3). But it can be easily checked that V"'(M C , 0) > for any 
value of 0i lying within subinterval IV. So from Theorem 9, we can conclude that there does 
not exist any solitary wave solution at M = M c for any value of 0i lying within subinterval 
IV. As the system supports coexistence of both NPSWs and PPSWs in subintervals II and 
III, i.e., for Pi lying within the interval 0i c < 0i < P 3c , from Theorem 10, we can conclude 
that there must exist solitary wave or double layer solution at M = M c when j3 lc < 0\< 3c 
and the nature of the solitary wave or double layer solution depends on the sign of V"'(M C , 0). 

In FIG. El M c , M D and V"'(M C , 0) have been plotted against fa for le < X < 3c . From 
this figure we see that the double layer solution is not possible for 0\ c < 0i < 03c, as the 
curve M = Md does not tend to intersect the curve M = M c for any point of 0\ lying within 
the interval 0\ c < 0\ < 0% c . So, for any 0\ lying within the interval 0\ c < 0i < 03 C , we get 
either a NPSW or a PPSW at M = M c . Again from FIG. El we see that the curve M = M D 
tends to intersect the curve M = M c at 0\ = 0i c , i.e., there always exists a PPDL solution in 
any right neighborhood of the point M c for some 0\ > 0\ c . Therefore, at 0\ = 0i c , we have 
a PPDL solution at M = M c . The existence of solitary wave solution at M = M c for any 
0i lying within the interval 0\ c < 0\ < 0^ c and the existence of PPDL solution at M = M c 
when 0i = 0ic can easily be verified through FIG. El Specifically, we have V"'(M c ,0) < 
for any 0\ lying within the interval 0\ c < 0\ < C and V'"(M C , 0) > for any 0i lying within 
the interval C < Pi < Ps c . Consequently, we have PPSW at M = M c for any 0i lying 
within the interval 0\ c < 0\< p c , whereas we get NPSW at M = M c for any 0i lying within 
the interval C < Pi < /3 3c and at 0i = C , there does not exist any solitary structure at 
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M = M C . 

To be more specific, we draw FIG. This figure shows that the amplitude of the PPSW 
increases with decreasing M c and ultimately this sequence of PPSWs end with a PPDL 
solution at M — M c when (3\ assumes the value f3\ c . However, we have seen in the literature 
that when all the parameters involved in the system assume fixed values in their respective 
physically admissible range, the amplitude of solitary wave increases with increasing M 
and these solitary waves end with a double layer of same polarity, if exists. Here it is 
important to note that M is not a function of the parameters involved in the system but is 
restricted by the inequality M c < M < Mp, where M = Mp corresponds to a double layer 
solution. So we cannot compare this case with the case of M c , since M c is a function of the 
parameters involved in the system and consequently, monotonicity of M c entirely depends 
on a parameter when the other parameters assume fixed values in their respective physically 
admissible range. But the solitons and double layer are not two distinct nonlinear structures 
even when M = M c . Actually, double layer solution, if exists, must be the limiting structure 
of at least one sequence of solitons of same polarity even when M = M c . More specifically, 
existence of double layer solution implies that there must exists a sequence of solitary waves 
of same polarity having monotonically increasing amplitude converging to the double layer 
solution, i.e., the amplitude of the double layer solution acts as an exact upper bound or 
Least Upper Bound (lub) of the amplitudes of the sequence of solitary waves. Therefore, if 
the double layer solution exists at M = M c then this double layer solution is also a limiting 
structure of a sequence of solitary waves of same polarity. So, double layer solution can be 
regarded as a source or sink of solitary waves of same polarity for increasing Mach number. 
Particularly, for the problem of Das et a/.—, the double layer solution at M = M c when 
01 = Pic can be regarded as a source for PPSWs for increasing M c along the curve M = M c 
whereas the same double layer solution at M = M c when j3i = 0\ c can be regarded as a sink 
of PPSWs for decreasing M c along the curve M = M c . However, for any other double layer 
solution, i.e., for the double layer solution at M = M D (> M c ) can be regarded as sink of 
PPSWs for increasing M restricted by the inequality: M c < M < Mp. 

On the other hand, if we take the interval (3 C < Pi < and move along the curve 
M = M c , then we have seen from FIG. H] that the amplitude of NPSW increases with 
increasing M c and it attains its maximum value at P\ = /?3 C . This is, of course, expected 
since /3 3c is the point of intersection of M c and M max and for X > /3 3c there does not exist 
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any solitary structure at M = M c . For this interval, i.e., ft < ft < ft c , the solitary wave 
solution at M — M c is similar to the solitary wave solution for M > M c , because for this 
interval the amplitude of NPSW increases with increasing M c . 

Moreover, it can be easily checked that the solitary structures along the curve M = M c 
is same as the solitary structures for M > M c , if one can move the solution space through 
the family of curves parallel to the curve M = M c . To be more specific, solution space for 
the present system with respect to ft has been presented in FIG. 0(a), in which the curve 
M = M c is omitted from the solution space as presented in FIG. 3 of Ref [26]. Now consider 
the family of curves parallel to M — M c . For instance, consider one such parallel curve 
for M = M c + 0.01 as shown in FIG. [5](a). In this figure, ft is the value of ft where the 
curve M = Mr, intersects the curve M = M c + 0.01, whereas ft is the value of ft where the 
curve M = M max intersects the curve M = M c + 0.01. FIG. 0(a) can be interpreted in the 
same way of FIG. 3 of Ref [26] with ft c replaced by ft and ft c replaced by ft. However, in 
FIG. [51(a) the solitary structures exist along the curve M = M c + 0.01, specifically, (i) only 
NPSW exists for < ft < /3 P , (ii) both NPSW and PPSW coexist for /3 P < ft < ft,, (iii) 
only PPSW exists whenever ft < ft < (3 M and at ft = ft, a NPSW coexists with a PPDL. 
The variation of amplitude of these solitary waves (exist along the curve M = M c + 0.01) 
have been shown in FIG.G^b). This figure shows that the amplitude of NPSWs increase with 
ft having maximum amplitude at ft = ft. Again, the amplitude of PPSWs decreases with 
increasing ft for ft > ft and the amplitudes of PPSWs are bounded by the amplitude of 
PPDL at ft = ft. Moreover, the amplitudes of the PPDL solutions increases with increasing 
Mo along the curve M = Mq. FIG. 0(b) also shows that along the curve M = M c + 0.01, 
the amplitude of PPSW increases with decreasing M along the curve M = M c + 0.01 and 
ultimately, these PPSWs end with a PPDL at ft = ft. This fact is similar for the case, 
when we move along the curve M = M c . So, if we move the solution space along the 
family of curves parallel to the curve M = M c , then one cannot distinguish between the 
solitary structures, exist along the curve M = M c and any curve parallel to M — M c . 
For ad = 0.0001, a = 0.95 and fi = 0.01, the values of ft and ft lie in the neighborhood 
of ft = 0.44475 and ft = 0.5809, respectively. To form a PPDL, the positive energy must 
dominate the negative energy in some region of the system and the PPDL will occur at some 
point where the potential drop is maximum. From FIG. E£b), it is clear that the positive 
potential dominates the negative potential in a right neighborhood of ft = ft and again the 
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potential drop is maximum thereat. This figure shows that physical reasons also suggest the 
existence of PPDL solution at fix = fi p , which is already confirmed in FIG. 0(a). This is the 
reason for the coexistence of NPSW and PPDL at fix = fi p . Actually, the same phenomena 
happens for the occurrence of PPDL solution at M = M c , although there does not exist any 
NPSW along the curve M = M c for < fix < fixc- We'll see later, that weakly nonlinear 
negative potential solitary structure exist along the curve M = M c for < fix < fixe and also 
the amplitude of the NPSW tends to zero as M — > M c , whereas PPSW of finite amplitude 
exist along the curve M = M c in the right neighborhood of fix = fixe and consequently, the 
potential drop is maximum at fix = fix c - This is the reason behind the occurrence of PPDL 
at M — M c when fix = fixc- This can also be interpreted as follows: Since there is no NPSW 
solution along the curve M = M c for < fix < fixc-, one can simply think the occurrence 
of NPSW of amplitude zero along the curve M = M c for < fix < fixc and in this case 
also the potential drop is maximum at fix = fixc and consequently, we have a PPDL solution 
at M = M c . So the solitary structures at M = M c is not localized. Actually, the solitary 
structure along the curve M = M c is same as the solitary structures for M > M c , if we 
move the solution space through the family of curves parallel to the curve M = M c . 

Again, at fix = fi ci V"'(M c ,0) = and V""(M c ,0) > 0, and consequently, = is the 
position of stable equilibrium. Therefore, it is not possible to get any solitary structure 
of the energy integral (20) of Das et al.,— at M — M c when fix = fi c . From FIG. El we 
see that along the curve M = M c the amplitude of PPSW decreases with increasing fix for 
file < fix < fic and ultimately it collapses at fix — fic- On the other hand, from FIG. HI we 
see that along the curve M = M c the amplitude of NPSW decreases with decreasing fix for 
fic < fix < fizc and ultimately it also collapses at fix = fi c . So, the figures FIG. [3] and FIG. 
H] are consistent with the physical interpretation for nonexistence of solitary structure at 
M = M c when fix = fi c . Again, PPSWs and NPSWs at M = M c behave KdV like solitons 
in the left and right neighbourhood of the point fix = fi c , respectively. 

Again, the existence of solitary structures at M = M c indicates the applicability of small 
amplitude theory to study weakly nonlinear solitary structures. To be more specific, if we 
replace V(<p) by V^{(f) to study weakly nonlinear solitary structures on the basis of the 
assumption \<p\ 1, the energy integral (20) of Das et ai.— assumes the following form: 




(23) 
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where V 3 ((p) is given by 

V 3 (<P) = X -A^ + X -B<\>\ (24) 

A - 1 - 1 ~ h + a/i (25) 

A ~M 2 -3a d 1-fM ' (25) 



B 



a 2 /j 3(M 2 + <r d ) 



2 [ 1 — A* (M 2 -3o- d ) 3 J' (26) 
We do not consider the case 5 = 0, since the expression for d 2 '(f) / 'd^ 2 , in this case, does not 
contain any nonlinear term and therefore no solitary structures exist. For B ^ 0, we can 
write V 3 ((p) as 

W) = f 2 (0 - 0m) = -|V (0 m - 0), (27) 

where 



3A 
2B' 



(28) 



Now applying the conditions for existence of the solitary structures of the energy integral 
([23]) as discussed in first para of Sec. II, we have found that weakly nonlinear DAS wave 
exists if and only if A < and B ^ 0. The condition A < gives 



M>M C = J3a d + — ^— . (29) 

V 1 + afjL — pi 

So, M c is the lower bound of M from which solitary wave solution of (123]) starts to exist. The 
condition B ^ determines the nature of the solitary wave solution of ([23]). In particular, 
the solitary wave solution of ( 12 3 j) describes either a PPSW or a NPSW according to .£> > 
or 5 < 0. The amplitude of PPSW or NPSW is given by \<p m \ = -(3A)/(2\B\) ->■ 
as M — )■ M c + 0. Furthermore, solution of (]23l is unable to describe the existence of 
both PPSW and NPSW simultaneously, and consequently, V 3 (4>) is representative of V(4>) 
to study weakly nonlinear solitary structures only when Pi lies either within the interval 
< 0i < I3i c or within the interval f3 3c < fii < (3 M since the amplitude of the solitary 
wave solution of f[2"3"]) tends to zero as M — > M c + 0. In FIG. El B is plotted against 0i for 
M = M c + 0.00001. Although, from FIG. |6](a), we see that B is equal to zero at 0i = (3 lc , 
one can easily verify that B becomes zero in the neighborhood of /?i = C when M is close 
to M c . Therefore FIG. |6]^a) and FIG. [6]^b) actually shows the existence of weakly nonlinear 
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NPSW or PPSW according to (3± lies either within the interval < ft < ft c or within the 
interval /?3 C < ft < ft./. It is also important to note that for < ft < ft c or ft c < ft < (3m, 
there does not exist any solitary structure at M = M c . 

To recover the coexistence of small amplitude NPSWs and PPSWs along with PPDLs 
when ft lies within the interval ft c < ft < ft c , we consider the following modified energy 
integral with n = 4. 

i ' " ^ 2 



Here V^(cf)) is given by 

V 4 (<P) = l -A<f + l -B<p* + hj<p\ (31) 



where 



3(5M 4 + 30M 2 a d + 9a 2 d ) 1 + 3ft + a 3 /i 



(32) 



(M 2 - 3<7 d ) 5 1 - /i 

A and .B are given by Eqs. (|25|) and ( 126|) . respectively. For C = 0, V^(0) is same as V 3 (0), 
so, we consider the case for C^O. For C ^ 0, we can write ^(0) as 
C 

V 4 ((j)) = — 2 (0 - 0ml) (0 - 0m2), (33) 

where 



^ml 



^m2 



(34) 



25 + - 18AC 



(35) 



2B - -18AC 
From ( )34|) and ( |35l) . we note the following facts. 

(i) 4> m i and m2 both are real and distinct if 4B 2 — 18AC > 0. 

(ii) m i0 m2 = 2A/C. 

(iii) (<p - m i)(0 - 0m2) < if and only if min{0 m i, 4> m2 } < < max{0 m i, m2 }. 

(iv) ml ->■ and m2 -> as M -> M c + 0. 

The necessary conditions for the existence of solitary structures of the energy integral ( |30l) 
are V 4 (0) = V^(0) and ^"(0) < 0. The conditions V A {0) = and V^(0) = are trivially 
satisfied whereas the condition V^"(0) < gives A < 0. Here also the condition A < gives 
M > M c . Therefore, as A < 0, from (ii), we see that m i0 m2 < if C > and in this case, 
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Vi(<j)) < for all min{0 m i, 4> m 2\ < 4> < as well as Vi(</>) < for all < < max{0 m i, TO 2} 
and consequently, both NPSW and PPSW exist simultaneously. On the other hand, if 
C < 0, one cannot get coexistence of weakly nonlinear NPSWs and PPSWs, since for 
C < 0, both <p m i and m2 are either positive or negative. Again as <f> m i — > and m2 _ > 
as M — > M c + 0, V±(<f>) effectively define the coexistence of both NPSWs and PPSWs when 
f3i lies within the interval (3\ c < fix < 3c provided that in this sub interval of fli, C > and 
AB 2 - 18AC > 0. In FIG. [3(a), C and AB 2 - 18AC are plotted against ft (f3 lc < ft < f3 3c ) 
for M = M c + 0.00001. A quick look of this figure may suggest that AB 2 - 18AC = 
at some point of ft. To remove this illusion, we draw AB 2 — 18 AC against (3i (in a small 
neighborhood of f3 c ) for three different values of M in FIG. [7(b). This figure clearly shows 
that AB 2 — 18AC — at fi\ — (3 C only when M = M c . However, we have already found 
that at (3\ = (3 C , <p = is a position of stable equilibrium, and consequently there does 
not exist any solitary structure of the energy integral at M = M c . Therefore, for M > M c , 
4B 2 — 18AC > and consequently, these two figures actually show the coexistence of weakly 
nonlinear NPSWs and PPSWs when 0% lies within the interval (3\ c < fl\ < 0.47(< Ps c ) . 
So, V^{(j)) is unable to produce coexistence of both NPSWs and PPSWs for all 0i within 
Pic < /^i < fizc, even when we take values of M smaller than M = M c + 0.00001. Again, in 
this interval of 0i, i.e., (3 lc < /?i < 0.47, the amplitude of both NPSWs and PPSWs tends 
to zero as M — > M c + 0. But for /3 lc < /?i < (3 C , we get PPSWs of finite amplitude at 
M = M c and for (3 C < 0x < 0.47, we get NPSWs of finite amplitude at M = M c . However, 
both PPSWs and NPSWs collapse at fli = (3 C , i.e., except for a very small neighborhood of 
(3i = (3 C , PPSWs having finite amplitude exist at M = M c and NPSWs of finite amplitude 
exist at M = M c . So, it is not possible to apply small amplitude method to study the weakly 
nonlinear solitary structures in a parameter regime where finite amplitude solitary structure 
occurs at M = M c , because in any small amplitude method, amplitude of solitary wave 
always tends to zero as M — > M c + which contradicts the occurrence of finite amplitude 
solitary wave at M = M c . There may be another possibility that V^cf)) is not the perfect 
match of V(4>) to study weakly nonlinear solitary structures in (3i c < (3\ < (3 3c - In any case, 
one can find inconsistency to study weakly nonlinear solitary structures for /3 lc < (3% < (5 3c . 
Therefore, our conclusion is that if there exists a region in the parameter space where finite 
amplitude solitary structures exist at M = M c , it is not possible to study weakly nonlinear 
solitary structures by means of any small amplitude method in that portion of the parameter 
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space. 



VI. CONCLUSION 

The paper deals with the existence of solitary wave and double layer at M — M c . The 
work is relevant to some recent numerical observations.— - — However, the existence of soli- 
tary wave at M — M c is limited by means of these numerical observations only. In this 
paper, our motivation is to find the clarifications behind these numerical observations. The 
analytical theory for the existence of solitary wave and double layer at M = M c has been 
provided for the first time in the literature. Moreover, this analytical theory is purely based 
on the physical interpretation of energy integral as discussed in Sec. II. Without having the 
prior knowledge of physics of energy integral, it is not possible to give the proper mathemat- 
ical clarification of the above mentioned numerical observations. Therefore, the theorems 
presented in this manuscript are based on the physical consequences of energy integral. 

On investigating DAS waves in SecJVl we have found the existence of double layer solution 
at M = M c . We have also discussed in our theoretical section, the conditions under which 
double layer solution occurs at M — M c . This observation is not specific for this problem 
only. Actually, in any plasma environments, if the curve M = tends to intersect the 
curve M = M c at some point of the compositional parameter space or solution space, then 
we can always get a double layer solution at M — M c , where M = Md{> M c ) corresponds 
to a double layer solution. We have mentioned this result without proof. The existence of 
double layer solution at M = M c has been observed for the first time in literature and thus, 
this is quiet a new result. 

From the investigation of solitary structures at M = M c of dust acoustic waves in SecJVl 
we find that the solitons and double layer are not two distinct nonlinear structures even 
when M = M c . Actually, double layer solution, if exists, must be the limiting structure of 
at least one sequence of solitons of same polarity even when M = M c . The following simple 
theorem of real analysis gives the actual logic behind the numerical observation that solitary 
waves end with a double layer, if exists. 

If {x n } is a monotonically increasing sequence of real numbers, then {x n } converges to its 

lub. Furthermore, if the sequence is bounded above, then lub{x n } is finite. 

This theorem gives the actual logic behind the numerical observation of Baboolal et ai.- 
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Actually, Baboolal et al— numerically observed that the double layer solution is the ultimate 
position of solitons, i.e., solitary waves end with a double layer, if exists. This numerical 
observation of Baboolal et al. 7 has been referred by several authors in different plasma 
environments. If we explain any fact with proper justification by means of mathematics, we 
can apply it in any specific problem. In the present paper, with the help of simple theorems, 
we have extended the physical interpretation for the existence of the solitary structures of 
the well-known energy integral even when the Mach number M assumes its smallest possible 
value M c . 
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(a) 




FIG. 1. Figure (a) corresponds to the subcase 1 of case 2 of Theorem 6. Figure (b) corresponds to 
the subcase 2 of case 2 of Theorem 6. Figure (c) corresponds to the subcase 3 of case 2 of Theorem 
6. 



a = 0.95 , n = 0.01 , a d = 0.0001 




FIG. 2. M c , Mb and V"'(M c ,0) are plotted against (3\ for values of Pi lies within the intervals 
Pic < Pi < /53c- At Pi = Pic, the curve M = Mp tends to intersect the curve M = M c . 
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FIG. 3. F(0) is plotted against for different values of /3\ lies within the intervals f3\ c < (3i < (3 C - 




FIG. 4. V(0) is plotted against <p for different values of j3\ lies within the intervals /3 C < /3i < /?3 C . 
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FIG. 5. (a)solution space for the present system with respect to ft. Here the curve M = M c is 
omitted from the solution space as presented in FIG. 3 of Ref [26]. (b) Variation in amplitude of 
NPSWs and PPSWs along the curve M = M C + 0.01. 




FIG. 6. B is plotted against ft for (a) < ft < ft c , (b) ft c < ft < /3 M - 
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FIG. 7. (a) C and 4B 2 - 18AC is plotted against ft for f3 lc < fa < fa c . (b) 4B 2 - 18AC is plotted 
against fa in a small neighborhood of fa. 
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